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1. Introduction and the main result 


In this paper, we consider the following density-dependent liquid crystal flows in 
' pt -h div(pw) = 0, 

{pu)t + div{pu (8) m) — pAu + VP = —V • (Vd 0 Vd), 

< dfd — 1/Ad-h M • Vd = I Vd|^d, |d| = 1, (1.1) 

divM = 0, 

^ p\t =0 = Pq{x), n\t=Q = uo{x), d|f=o = do(x), |do(x)|=l in R^ 

where p is the density and ii is the velocity field, d stand for fhe unif vecfor field which represenfs 
the macroscopic orientations and p an v are two positive viscosity constants. The term Vd O Vd in 
the equation of conservation of momentum denotes the 3x3 matrix whose (i,/)-th entry is given 


by did ■ djd (1 < z,/ < 3). 

When d is a given constant unit vector, dl.lll is reduced to the well-known inhomogeneous in¬ 
compressible Navier-Stokes system, which has been studied by many researchers. When the vis- 

n 

cous coefficient equals some positive constant, Ladyzenskaja and Solonnikov 112211 first established 
the unique resolvability in a bounded domain with homogeneous Dirichlet boundary condition for 
u; similar result was obtained by Danchin 191] in R" with initial data in the almost critical Sobolev 
spaces; Simon 13411 proved the global existence of weak solutions. In general, the global existence 
of weak solutions wifh finite energy with variable viscosity was proved by Lions in (see also 
the references therein, and the monograph fj]). Yet the regularity and uniqueness of such weak 
solutions is a big open question in fhe field of mathematical fluid mechanics, even in two space 
dimensions when the viscosity depends on the density. Except under the assumptions: 


Hipo) 


po e L“(T2), inf „ 

Oo c 


— 1||l”(t 2) < L and uq G H^(T^), 


Desjardins 11311 proved that Lions weak solution {p,u) satisfies u G L“((0,T);H^(T^)) and p G 
L“((0, T) X T^) for any T < oo. Moreover, with additional assumption on the initial density, he 


could also prove that ii G L^((0, t); H^(T^)) for some short time T. Gui and Zhang 11611 proved 


the global wellposedness of inhomogeneous Navier-Stokes system with initial data satisfying ||po ~ 
l||^s+i being sufficiently small and mq G H®(R^) n H^'^(R^) for some s > 2 and 0 < e < 1. However, 
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the exact size of |||0o ~ not given in [16]. Very recently, Danchin and Mucha il2ll proved 

that: given initial density po G L“(Q) with a positive lower bound and initial velocity Uq G 
for some bounded smooth domain of IR”, the Navier-Stokes system with constant viscosity has a 
unique local solution. Furthermore, with the initial density being close enough to some positive 
constant, for any initial velocity in two space dimensions, and sufficiently small velocity in three 
space dimensions, they also proved its global wellposedness. When the density p is away from 
zero. In llO, Abidi proved in general space dimension n that: if 1 < p < 2n, 0 < p < }i{p), given 


flo G and wq G ^(R”) inhomogeneous Navier-Stokes system has a global solution 


-l+n/p. 


provided that HflolLn/p + ll^olL-i+n/p < Co for some sufficiently small Cq. Moreover, this solution is 
unique if 1 < p < n. This result generalized the corresponding results in [8,9] and was improved by 
Abidi and Paicu in |3] with Uq G and Uq G (M.") for p,q satisfyin g s ome technical 

assumptions. Abidi, Gui and Zhang removed the smallness condition for uq in H]. Notice that 
the main feature of the density space is to be a multiplier on the velocity space and this allows to 


define the nonlinear terms in the system. Recently, Danchin and Mucha lllll proved a more general 


wellposedness result of this system with p(p) = p > 0 by considering very rough densities in some 
multiplier spaces on the Besov spaces for 1 < p < 2n, which in particular completes 


the uniqueness 


motivated by llTl. 


3 result in [10 for p G {n,2n) in the case when p(p) = p > 0. On the other hand. 


32l. 1381] concerning the global wellposedness of 3-D incompressible anisotropic 


Navier-Stokes system with the third component of the initial velocity field being large, Paicu and 


Zhang 13311 proved that: given Aq G B”^‘^(R^) and Mq = {uq,Uq) G B^^ ^(R^) for 1 < q < p < 6 


1”^(R3) and uq = G 

and 1/q — 1/p < 1/3, with p(p) = p > 0 has a unique global solution as long as 


p||flo||g3/? + I|woIIb-]+3/p) exp I Co 11 Mol 1 


for some sufficiently small Cq. We emphasize that the proof in [3311 used in a fundamental way the 


algebraical structure of the momentum equation, namely div u = 0. 

When the density p is a constant, [CT reduce to be the classical nematic liquid crystal system 


which was first first introduced by Lin [2611 as a simplification of the Ericksen-Leslie model il4l. 


23D. 


As it is difficult to deal with the high nonlinear term | in the third equation of dl.ll) . Lin and 


Liu |l26t 


proposed to study an approximate model of liquid crystals system by Ginzburg-Landau 
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function, i.e., \Vd\^d is replaced by —f{d) = —VF{d) = — V-^(|d|^ — 1)^ with e > 0, and the third 
equation of dl.ll) is replaced by dtd + (m • V)d = }i{Ad — f{d)). In this situation, Lin and Liu 12611 


established the existence of global weak solution in dimensions 2 and 3 and assumptions that uq G 
and do S . Existence and uniqueness of global classical solutions were also obtained by them 
in dimension2 provided Uq G and do G and provided the viscosity v is large in dimensions. 


The partial regularity of suitable weak solutions was studied in [271]. However, as pointed out in 
1^ . the vanishing limit of e ^ 0 is an open and challenging problem. For the full system , in two 


independent papers [2011 and [25(1 . the authors respectively established the global existence of Leray- 


Hopf type weak solutions to liquid crystals system in dimension2, and proved that the solutions are 
smooth away from at most finitely many singular times which is analogous to that for the heat 
flows of harmonic maps (see lO]). The global existence of weak solutions in dimension 3 is still now 
an open problem. For strong solutions, the global existence of strong solutions with small initial 
data and the local existence of strorw solutions with any initial data of system have been studied 
by Fi and Wang |24], Lin and Ding |29], Wang |35|, Hineman and Wang |19], and Hao and Liu 118|. 


Recently, the uniqueness of weak solutions to this system has been studied by Lin and Wang [2811 


Blow up criteria for local smooth solutions has been studied by [2ll. Bill and the references therein. 


Very recently. Wen and Ding [36 ] obtained the local strong solutions to the problem (11.11) with 


vacuum. Fan, Zhou and Nakamura [iSH extended it to be global in 2D in a bounded smooth domain. 


In this paper, we suppose that the initial density verifies inf po (^) > 0 and thus, by the maximum 
principle for the transport equation, we have inf p{t,x) > 0. We define a = 1/p — 1 and assume 
p = V = 1 for convenience which allows us to work with the following system: 

' dtu + u ■ Vfl = 0, 


dfU + u- VM + (l + fl)(VP — Am) = — (1 + fl) V • (Vd © Vd), 

< 3td — Ad + u ■ Vd = I Vd|^d, |d| = 1, (1.2) 

divM = 0, 

^a\t^o = ao{x), u\t^o = uo{x), d\t=o = do{x), |do(y)|=l in R^. 
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The main result of the present paper is stated in the following theorem: 


Theorem 1.1. Let uq G {uo,do) G ^ divwo = 0/ 

1 + uq > b 


(1.3) 


for some positive constant b. Then there exists a positive time T such that (lf.2D has a unique local-in-time 
solution {a, u, d) with 


a G Q([0,T];B3/2(r3)^ u G Q([ 0 , T]; b1/2(]r3) ^ 
d G Cb([0,T];B3/2(R3) nLi([0,T];B2"/2(R3)). 

Moreover, There exist a positive constant C such that for any {ao, uq, do) verifying 


C 


|lflo|lB3/2(l + ||MoIIb1(2 + (I|MoIIb1/2 + II^o|Ib3/2)^) < 


(1.4) 


and 


C(||woIIb1/2 + IMo|Ib3/2) (1 + (II“oIIb1:;2 + II‘^o|Ib 3/2)^ + ll“ollgV2(l|wollBi:;2 + II^o|Ib3/2)^^^ ) < 1. (1.5) 
2,1 2,1 \ 2,1 2,1 1^21 2,1 2,1 

Then the system (ll.ip has a unique global solution {a, u, d) with 


fl G C([0, + 00 ); Bl'^) n L“((0, +oo); 

U G C([0, +00); n £00£1/2) p £5/2)^ 

d G C([0,+Oo);B3/2) p£oo^q^^^).£3/2) 

The remainder of this paper is organized as follows. In section 2, we shall first review some 
useful statements on functional spaces and basic analysis tools, and introduce several technical 
Lemmas. In Section 3, we give the estimates of the transport equation. In Section 4, we shall prove 
the local well-posedness of Theorem 1.1. In last section, we complete the proof of the global well- 
posedness. 

Notations : Let A, B be two operators, we denote [A; B] = AB — BA, the commutator between 
A and B. For a <b, we mean that there is a uniform constant C, which may be different on different 
lines, such that a < Cb. For X a Banach space and I an interval of R, we denote by C(f; X) the set 
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of continuous functions on I with values in X. For cj G [1, +oo], the notation U{1; X) stands for the 
set of measurable functions on I with values in X, such that t ||/(f) ||x belongs to For a 

vector V = {vi, V 2 ) G X, we mean that all the components Vi{i = 1, 2) of v belong to the space X. We 
always denote is a generic element of /^(Z) so that dy = 1. 

2. Preliminaries 

Let {x, (p) be two smooth radial functions, 0 < {x, (p) < 1, such that x is supported in the ball 
B = {^ e 1^1 < 4/3} and (p is supported in the ring C = {^ G R^,3/4 < |^| < 8/3}. Moreover, 
there holds 


Y2 ^ 0. 

jez 

Let h = and h = T~^Xf then we define the dyadic blocks as follows: 


V = <p(2 ’D)f = h{2Jy)f{x - y)dy, 

JR3 

Sjf = x{2-^D)f = [ h{2iy)f{x - y)dy. 

JR3 


By telescoping the series, we thus have the following Littlewood-Paley decomposition 


u= Yj^jU, yu E S'(R^)/P[1R.% 

where /^[IR^] is the set of polynomials (see|l3)- Moreover, the Littlewood-Paley decomposition sat¬ 
isfies the property of almost orthogonality: 

AjtAyu = 0 if |fc — ;| > 2 and Ajt(Sy_iMAyM) = 0 if |fc — j\ > 5. 

Now we recall the definition of nonhomogeneous Besov spaces. 

Definition 2.1. Let s < 3/p (respectively s G Rl, {r,\,p) G [1,-|-oo]^ and T G (0,-|-oo]. We define 
Lj(Bp,,(R^)) as the completion o/C([0, T];iS(R^)) by the norm 


L}(Bp 


%.) = L 2'” 

Laez 


'0 


j \ r/A' 


1/r 


< 00 , 


zvith the usual change ifr = 00 . For short, zvejust denote this space by 1^(6®,,). 
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Remark 2.2. It is easy to observe that for p, r, X, X\, A 2 G [1, + 00 ], 9 G [0,1], 0 < Si < $ 2 , we have the 
following interpolation inequality in the Chemin-Lerner space (seelsl]): 

with 1/A = 0/Ai + (1 — 0)/A2 and s = 0Si + (1 — 6 )s 2 . 

Let us emphasize that, according to the Minkowski inequality, we have 

ll“ll4(B“,,) — II^^IIl^(bs_,) if > L I|w||£A(bs^^) > ||w|Il^(bs_^)/ if a < r. 

The following Bernstein's Lemma will be repeatedly used throughout this paper. 

Lemma 2.3. Let B be a ball and C a ring o/R^. A constant C exists so that for any positive real number A, 
any non-negative integer k, any smooth homogeneous function a of degree m, and any couple of real numbers 
{a, b) with 1 < a <b, there hold 

Suppii CXB^ sup Wd^uW^i, < 

\oc\=k 

Suppii C AC => < sup ||9 “m||l“ < 

\a\=k 

Suppii C AC => \\cr{D)u\\i^b < 

In order to prove the global well-poseness of our main Theorem, we also need the following 
anisotropic Bernstein's Lemma . 

Lemma 2.4. (see^) Let Bh (resp. Bv) be a ball o/R^ (resp. Ri,) and Ch (resp. Cv)a ring o/R^ (resp. Rd). 
Let 1 < p 2 < Pi < 00 and 1 < q 2 < Then there hold: 

If the support off is included in 2fBh, then 


If the support off is included in 2 By, then 

If the support off is included in 2^Ck, then 






LiHLi^y 


sup ||a“,/|lBP 


|a|=N 


^hJ wqHiiy- 
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If the support off is included in l^C^, then 




In the sequel, we shall frequently use Bony's decomposition from in the homogeneous con¬ 


text: 


uv — tuV + %u + R{u,v) = tuV + TZ{u,v), 


( 2 . 1 ) 


where 


TuV = Sj-iuAjV, R{u,v) = 


and 


hjV= ^ Ay/z;, n{u,v) = Y,Sj+2oAj 

li-y'l<i yez 

In the sequel, we shall frequently use Bony's decomposition from in the homogeneous con¬ 


text: 


uv = tuV + R{u, v) = TuV + tyU + lZ{u, v), 


where 


tuV = Y^Sj-iuAjV, R{u,v) = AjuSj+ 2 '^, lZ{u,v) = Y, AjuAjV and AjV = Y 
jez jez jez 

As an application of the above basic facts on Littlewood-Paley theory, we present the following 
product laws in Besov spaces, which will be constantly used in the sequel. 

33/) Let Si < 3/2, S2 < 3/2 with Si -I-S2 > 0. Let a G B 2 ^{Rf),b G Then 


Lemma 2.5. (see/y, 
ab G 


||fl&|Lsi+S2-3/2 < C||fl||gSl ||fe||g=2 . 

^2 1 2,1 2,1 


Lemma 2.6. (seel3]) Let si < 3/2, S2 < 3/2 with Si -|- S2 > 0. Assume that f G B 2 \(K^) and g G 
^ 2^00 (^^)- Then there holds 

ll/S'llBn+^2-3/2 < C||/||g^l ||y||B^2 


^2,co 


^2,1 ^2,co 


Lemma 2.7. illd/ Let s G (—3/2,5/2]. There exists a sequence Cq G i^{^) such that ||ci||/i = 1 and 


constant C such that 


\\[v-V,Aq]u\\q^2 <CCq2 ‘^1|Vz7||g3/2||M||B; 

2,1 


( 2 . 2 ) 
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In the limit case s = —?>I2, we have 


sup2 -V,Kj\u\\i 2 < C||Vz;||g3/2||wL-3/2. 

;eZ 

To deal with the pressure term, we also need: 

Lemma 2.8. [2] (i) Let a G and VP G L^(R^), then there holds 

\\[a,Ac]VP\y < Cc,2-‘?||fl||g5/2||VP|L2. (2.3) 

(ii) Let a G H^(R^) and VP G then there holds 

||KA,]VP||,2 < Cc^2^/^a\\H4VP\\H-r. (2.4) 


The proof of the following two lemmas can be found in 1371] . here we also outline its proof for 
completeness. 

Lemma 2.9. For 2 < q,r < oo,v = G L“((0,+oo);B^^^) n L^((0,+oo);B^^^) with divi? = 0, 


there hold 


1^^ IIl?(L’( 4)) < Cd;2 llg2(g3/2) 11^^ 


(2.5) 


and 


\AjV^ 


'LpLliLD) 


^ IW II^Un5/2^ 11^ 


,h\\l/2-l/r 




( 2 . 6 ) 


Proof : Firstly, according to divi? = 0 and the following version of Gagliardo-Nirenberg inequal¬ 
ity: V/ G 22 (R) 


(2.7) 


we have 

WAjV^Wll < C||Ayr;3||[^2+l/r||9^^,^3||l/2-l/r 

< C|lAyr;3||i/2+i/>-||A.(div;,r;'')||[^2-iA < C2Ai/2-i/>-)||A^-r;3||i/2+i/r||^.y/!||i/2-i/r_ (2.8) 

Thus, by Lemma l231 we have 

< (^2A3/2-2/<?-l/d||A.„3||l/2+l/r \\k hul/2-l/r 

< Cd-2“A2/i?+l/d||7;3||l/2+l/n| ^||l/2-l/r q n\ 

lit; ||^(g3/2)ll^ llt2(g3/2)- (^-V) 
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Similarly, we have (I2.6D . From the above Lemmas, we can obtain the following important Lemma. 


Lemma 2.10. Let u G L’^{bYi) H LJ(B 2 i^) div u = 0, one has 


2^||Ay(wV)||Li(L2) + ||Ay(w3div;,M'A||Lj(L2) 


<cdi 2 iiTr.,t,i/z 




“*■(11^^ llLi(Bff) IIl“(b1_(^)) (II“ IIlKbI^) 11*^ llL“(B^f)) 

=Cd^2-t'^F{i? ,u^). 


y3/2 


( 2 . 10 ) 


Proof : Firstly, thanks to Bony's decomposition, we have 


+ U^diVhU^ =7’u 3M^ + + r„3diV;,i/^ + 

+ R{u^,u^) + R(m^, div/jW^A- 


( 2 . 11 ) 


Applying Lemma l23l and (12.51) with q = r = oo gives rise to 


2lAy(T„3M')||,j(,.) + ||Ay(r,3div,M')||,j(,3) 

< C ^ (2lSy,_iM3||^2(^^)||Ay,M'^||^2(L2) + || ||^2 (l~) || AyAdiV/^w'^ IIl?(L2)) 

\f-i\<5 

<C ^ (^2AAy//lA|A2(£c.)||Ay/M''|A2(£2) + ||Ay//M^|A2(^oo)||Ay/(diV;,i/^)|A2(^2)] 

||M IIi2(b3/2)I|W 1^2(83/2) 


A/2 


( 2 . 12 ) 


Similarly, we have 


H E (2^1V(di'^^“^)llL?(L3)ll4'M^llL3(LAL~)) 


l/'-; 1 < 5;"<;'-2 


< Cdj 2 -i^^\\ 


.3||l/2 


lL?(BiA)'l“ "l?(B 2T) 


,hi\3/2 


< 


Cdj2 ^^^(||w^lAi(B5/2y + ll“^llL“(BAi^)^^^^^ll^^^'llLHBfA) “*■ ll“^llL“(BAih'^^^^‘ 


(2.13) 
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And while (I2.6I) applied with q = 2,r = oo gives 


l/'-; 1<5 
l/'-; 1 < 5;"<;'-2 

-Cdjl ^b 5/2) 11 11 ]. 1 ^g5/2) 11 w'i Ir ) 


(2.14) 


and 


2^||AyR(w^,M'')||gi(g2) <C ^ f2^’||Ay/l/^||£c<,(g»(g2))||Ay/l/^||gi(g2(g»)) 

j'>hNo 


<C ^ f2^+^ ||Ay/M ||L»(L2)||Ay/W ||gi(g2(g»)) 

7'>;-No 


-Cdjl ^^^\W^\\]i ^b5/2) 11 11 ]. 1 ^b5/2) 11 m '' 11 (b1_/2 


)• 


(2.15) 


In what follows, we give the estimate of the remaining term R(m^, div/jM^'). Applying (12.51) with 
q = 2, r = 00 gives 


l|A;R(MAdiV;,M^)||gl(g2) 

<C ^ (2^'||Ay,(div,,W^)Ay,u3||gi( l(g2^) 

<C ^ (2^||Ay,(diV;,w'')||g2(g2)||Ay,u3||g2(g2(L^)) 

;'>;-No 

< IIm^ II-2?■3/2^ IN^ II-2^3/2n 

- I II llL?(B3f)ll llL2(B3f) 

< Cdjl {\\u |Ib1(b5/2^ + ||m llL~(B 2 f)) ^ ^11^^ I 


Li(Bff) + ll^'^llLr(Bif))^^^ 


(2.16) 


Combining with the estimates (I2.13ll - (I2.16I) . we can finally get 


2^’||Ay(M3M'A||Bi(B2) + ||Ay(M3div,,M^)||gi(g2) 


<Cdj2 M|l/ ||gi(B5/2)ld^ llLi(B5(h^' lllriBzf) 


+ 


l“^llL)(B®f) “*■ ll“^llL”(Bif)) ^ (ll“ llLi(Bfj^) + 11“ IIl“(b1(^)) 


\ 3/2 


(2.17) 


The proof of Lemmais complete. □ 


3. Estimates of the transport equation 


Let us first recall standard estimates in Besov spaces for the following linear transport equation: 
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(3.1) 


dtf + v-Vf = g, 

/(^.O) =/o. 

Proposition 3.1. (see o) Let v be a divergence free vector field with Vv & L^([0, T]; Lor 

s G {—5/2,5/2],givenfo G G L^([0, T];B|^(R^)), tBen the system (13. IP has a unique solution 

f G C([0,T];B|^(R^)) . Moreover, there exists a constant C such that for \/1 G [0,7] 




c 




(3.2) 


Proposition 3.2. [2] Let m G ^ B 2 ^^(R^),Vm G Li^(B 2 ^^) with divw = 0, and a G 

C([0, T];B 2 (^(R^) such that {a,u) solves 


dtu + u ■ Vfl = 0, 
fl(x,0) = flo- 


(3.3) 


Then there hold for \/t <T 


< ||flo|lB3/2e 


.CU{t) 


\\a - Smfl||£~(B3/2^ < Y2 + ll«o|lB3/2(e' 

' q>m 


CU{t) _ 


1). 


(3.4) 


(3.5) 


with U{t) = ||VM||j^i(g3/2y 

Finally, we give estimates for the following constant coefficient parabolic system: 

dtu — aAu = f, 
u{x,0) = Mq. 


(3.6) 


Proposition 3.3. (see nl(jl £(c) Assume that cr > 0, then there exists a universal constant k such that for all 
s G R and T G R+, 

II“IIl~(b|_j) ^ I|wo|Ib|_j + II/IIli,(b|j)' 




— e 


-kcTA'IT 


)(II^<?“o||l2 + \\^qf\\L\{L2))- 


4. Local wellposedness of Theorem 1.1 

4.1. Existence of Theorem 1.1 

The proof of existence of a solution is performed in a standard manner. We begin by solving an 
approximate problem and we prove that the solutions are uniformly bounded. The last step consists 
in studying the convergence to a solution of the initial equation by a compactness argument. 
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Step 1: Construction of smooth approximate solutions. 


We first smooth out the initial data. For m G N, let flg = Wg = S„uo — S_„mo, 


dg = Sndo — S-ndo then (Aq,U g^dg) G H“(R^). Now, owing to |36], we deduce that system dl.lll 


with the initial data (aq, Wq, dg) admits a unique local in time solution (a", u", d'^VP") verifying 


a" G C([0,r];H®+i(]R3))^ 
u” G C([0,r];H^+i(K3))nLi„,([0,T”];H^+2(jj^3))^ 
d” G C([0,r];H®+2(]R3))nL?„^([0,r];H®+3(K3)) 

and 

VP” G L1([0,T”];H"(R3)) withs > 1/2. 

Step 2: Uniform estimates to the approximate solutions. 

Our goal is to prove that there exists a positive time 0 < T < inf T”(n G N), such that 
(a", u”, d”, VP”) is uniformly bounded in the space 


= ^ 10 




^2,1 


n X (^ 2 ^) n X 


1/2', 
2,1 / 


Let u’^ = u’[ + u",d" = d"+ d" which w” = e^^UQ,d’l = e^^dg, then {a",u",d",'VP") solves 

' dta" + {ul + u") ■ Va” = 0, 

dtu'^ + ul ■ Vm” + (1 + a”) (VP” - All”) - Mn, 

< dtd" - Ad'” + ■ Vd'” = G„, (4.1) 

divw” = 0, 

, (A”,w”,d'”)|f=o = (flO/0,0), 

with 


Mn=- ul ■ Vul - m” • Vul - t7” • Vt7” + a^Aul 

+ (1 + a"){Vdl ■ Adi + '^dl ■ Ad" + Vd'” ■ Ad'” + Vd'” • Adi) 


G„=- ul ■ Vdl - u" • Vdl - u" ■ Vd'” + iVd^j^d)) 

+ \Vdl\^d" + IVd'”|2d^ + IVd'”|2d'” + 2Vdl ■ Vd"dl + 2Vd^ ■ Vd'”d'”. 
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It is easy to observe from Proposition |33] that 


and 


Let 


^‘'L\\lf={U+;Bl{^) — ^ll“ollB2f' II^LllL”(R+;B®f) — ^II^o|Ib3/2 (4.2) 

IWIlL|(8>f) sex; 21«(1 -e-'«’')||A,u„||„, 

qez 

||dEllLj(B,^f) < C E - e-^^^"^)\\A,do\y. (4.3) 


Z”(f) =||m” 


Lri^Vi) 





+ ||VP"||,.,s.„, 


Firstly, we give the estimate of d. 

Applying to the third equation of (I4.1D , and then a standard commutator's process, gives 


dtA^d" + III • VAgd" + iZ" • VA,d" - AA^d" 

=[m 2 • V,Ag]d" + [iz" • v,Ag]d'^ + A^dv^l^o 
+ A, (I Vdl\^d ^) + A, (I Vd” \^dl) + A^{\ vr ) 

+ A^{2Vd’l ■ Vd"dl) + A,(2V< ■ Vd"d") + Ag(i/;( ■ Vd’l) + A,(iZ'’ • Vdl). (4.4) 

Taking the inner product of the above resulting equation with A^cf” and integrating with respect 
to t, we can finally get 


< E V,Add"lLi(L2) + E • V,AddlL|(L3) 

ijCZ ijCZ 


+ IW ■ v<ll,.,s», + llfl" ■ + |||v<|^d“||,.|s„, 

+ |||V<i"|V;||^|j,„, + |||Vi”|V”||,,,j„, + ||V< ■ + ||V< ■ ( 4 . 5 ) 


Applying Lemma 123112.71 yields that 


E2V2||[„«.v,Add'l,i(,2)<C 


“lIIb5/AM'1Ib3/2^t^ < C 


MEIlB5/2Z”(T)dT, 


(4.6) 


E 2 ^^/^\\[u" • V,Add"|Ll(, 2 ) <C [ ||ti'd|B 3 / 2 ||d'^||, 3 / 2 dT < C{Z"it)f 

qez. -^0 


(4.7) 
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|Ml- V<||^l(g3/2) <C^ l|WLllB3/2|MLllB5/2dT< Will 

II^^lIIbi(2||wEIIb5/2^'^ +II'^lIIb^/zIMlII 

<C(||Wo||b1/2 + ||4|Ib3/2)(||w"||b1(b5/2) + Wdl\\lll^Bl{^)) 


Lll]5V2l|MLllgV2lKll]5V2lMLllg7/2rfT 
^2,1 ^2,1 ^2,1 ^2,1 


i/idT 


(4.8) 


\u" -Vdl 


|Ib1(b 3/3) <C^ ||M-"||g3/2||d2||B5/2dT< l|M"l|-g2l|Ml|- 

<C(Z'^(0)2 + C||do|lB|/2||d2||Bj(B^7/2), 


]l5%\\dlW]l3%\KWg;%dT 

^2,1 ^2,1 ^2,1 


(4.9) 


\\\^‘>L\"‘>L\ln,im. <cf \\dl\\„^\dl\\l,„dT <cf \\dl\\l„\\dl\\^,ndT < C||rf„||j3,.ll<llt.,e;C)- 

t V 2,1 ^ JO 2,1 ^ 2,1 JO 2,1 2,1 ^ 2,1 2,1 ^ 

(4.10) 


|||V<|2d”||,j(g3/2) <cj^ ||(f'l53/2||<|||v2dT 




|||Vd'f<||,4(^3/2) <C ||<||B3/2||d'^|||3/.dT 


<C / ||<||B3/2||d'^||B3/2||d"||B7/2dT<C||do|lB3/3(^^ 


ll|Vrf”Prf1lB|(B3/3) <c/ l|d" 11,3/3 ||d-III, 


dr 


<c 


./' 


7)11|2 
B: 


3/2||rf1lB3/3rfT<C(Z'’(f))l 


' 2,1 


^2,1 


(4.11) 


(4.12) 


(4.13) 


WVdl-Vd’’dlW^,^^3n^<C / ||<||B3/2||J'lB5/2||<||B3/2dT 


^”'''{/\||d1l|02||<fj2||<irj2rfT 


'B 


<c 




33/2 


n7/ldT + C 


2,1 

<C(Z'’(f))'+C||do|||3/2||<||Bl(B,3/2 


I D/ 

'2,1 ^2,1 

II^^Lll63/2||df,|U7/2dT 


/o 


*LIIb|/2||<^LIIb3j 


livd; ■ vJV'iimss, <cjl ||rfj||,|,2||<i»|l8j,=||<i"|isWT 

SCj(Vzils£ll<ll-|ll‘f|l|£ll<i"ll-^ 

<Cj('||<i;|l8wll<i”Jl87«<iT + C Jl \\d"\\l,„m,,ndT 

<C(Z"(())‘‘ + ||<io|l8wl|<i;il8.,8;,. 


||7|2^T 

^2,1 


(4.14) 


(4.15) 
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Taking above estimates into (I4.5I) , we have 


<C(l+||rfo||B3,0(Z"(0)" + C(Z"(0)* + C^'(||<||j5,. + (l+||d„||j,«)IKIlB7,.)Z"(T)dT 

+ C fllwollgi/z + (1 + ||^^o|||3/2)(l + IMo|Ib 3/2)^ (II^LlIlKBlf) + 


In what follows, we mainly discuss the momentum equation. In the general case of large pertur¬ 
bations of a constant density state, solving the nonconstant coefficients elliptic equation may be a 
complicated problem. Here, we follow the Proposition 6 in to give the following key proposition: 

Proposition 4.1. Let a € L“(B2i^)(3 ^^)/wq S which satisfy 1 + a > b for some positive constant 

band divwo = 0. Let f belong to Lj{bI^^^) with divf G Lj{H^^),v G Lj(By^^),'Vv G LjlB^^-^) and 
{u, VP) G C([0, T]; n T); bI'^) x L]J{0, T); bI{^), which solve 


dtu -|- V ■ Vw — (1 -|- a){Au — VP) = /, 
divM = 0, 
u{x,0 ) = Uq. 


(4.17) 


Then there holds for Vf G (0, T] 


W^WipiBlf) + ll“llLi(Bff) + W'^^WlliBlf) 

^ C^||Mo||gi/2 ||/||gi(gi/2) +2^'"l|Sm«||L”(L2)||div/||^i(^_i^ 

+ l|M(^)llBi/2||t7(T)||g5/2dT +W(f)^ ||M(T)||gi/2dTj, 


(4.18) 


with 


^(0 ll'Smfl||g~(£2) (1 + ||®|Il~(l~) + II^IIl“(b 1/2P + 2 lkllL“(L°«)/ 

provided that 

\\a — Smfl||£~(g3/2^ < c 

for some sufficiently small positive constant c and some integer m G Z. 
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Applying Proposition l4.1l to the second equation of Il4.ll) . we have 


with 




+ 


l|M"(T')llBifl|WL(^)llB5/2^^ +||M"(T)||gi/2dTy 




M„(t) =-ul- Vul - a" ■ Vul - a" ■ Vu" + a^Aul 

+ (1 + a"){Vdl ■ Ad’l + Vdl ■ Ad” + Vd” • Ad'” + Vd'” • Adi), 


under the assumption 


|/7 SmU |||^oo^g3/2^ ^ 


Noticing that ||/^||h-i < C||/||gi /2 ||^||bo^ and 


we have 


and 


ll/'Vjsliy <c||/||y||Vg|L. <c||/||j.,.||g||'ft|k|l5«, 

2,1 £>2_i 02,1 

I|div(-M^ • Vul - w” ■ Vul - fl” • Vfl” + a’^Aul) ||^-i 

< ||,,n ||5/4 II n ||3/4 _i_|L7)i|| lli;’U|l/4 ||3/4 
~ INL II nl/2 II^L II d5/2 +11^ IIbH^ II^L II d1/2 Iu^L II d5/2 

^2,1 ^2,1 2,1 Djj D2,1 

2,1 0^1 0221 


n||l/4 ||,,«||3/4 _|_ ||,7»J||5/4 ||,7»J||3/4 
L II d1/2 II “L II d5/2 “T II “ IId 1/2||“ lln5/2 
^21 ^2^^ £>2,1 1 


||div((l + a''){Vdl ■ Adi + Vd^ • Ad'” + Vd'” • Ad'” + Vd'” ■ Ad';'))||^-i 


<(i + II«"IU3/0(K 


B3/2||«L|| A/2 + 11 «L 11 63/2 
021 2>2,i 2.1 


,,, ,.133/211^^ IId7/2 

2,1 Bji B 21 


2,1 021 021 021 021 


(4.19) 


(4.20) 


(4.21) 
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Combining with (I4.20D and (I4.21I) and applying the Young's inequality yields that 


2 "'||Smfl”||Lr(L2)||divM„||Li(H-i) 


+ Cf2®'"||S^fl"||^^^j^2)(l + ll'^”ll£~(B3/2p(l + ||Mo|||i/ 2+ ||do|||3/2) 

+ Cf2®'"||S^fl"||^^(j^2)(l + ||«”||£»(g3/2p(||wo|lBi/2 + ||do|lB3/2)(Z"(t))'^ 

Applying product laws in Besov spaces gives 


(4.22) 


^ hi ||bV2 hi IIb5/2 + ||m” llgi/2 ||w” ||g5/2 + ||w” llgi/2 ||Wl ||B5/2dT 

+ (1 + 11 fl" 11 gc. (b3/2) ) ^ (11 Ml 11 b|/2 + 11 Vdl 11 g3/211 11 gi/2 + 11 Vdl \ \ ^ 3/2 \ \ Ad" \ \ gi /2 

+ II V(P ||g3/2 II Ad” llgi /2 + II Vd” ||g3/2 II Ad” llgi/ 2 ^ dr 

~i) ll“”llBlf ll“”llBff + l|M”|lBV2||MEIlB5/2dT 

+ (1 + ||«”||g~(B|/2p hiyn + WdlWl^n + ||d”g||g3/2||d”||g3/2 + ||d”||5,,2dT 


< 


MlIIb 1 / 2 ||MlIIb 5/2 + ||w”||gl/ 2 ||t 7 ”||g 5/2 + ||M”|lBi/ 2 ||MLllB 5 / 2 dT 


^2,1 


^2,1 


^2,1 


2,1 


^2,1 


^2,1 


+ (1 + ||«”||g.(B3/2p [hiyn + ||d”g||g3/2||d^||g2/2 + ||d-”||g3/2||d”||g2/2jdT 
<(1+ ||fl”||gc.(g3/2))(||Mo|lBV2 + l|do|lB3/2)(||M^||gl(g5/2) + 11 d” 11 gl (g7/2) ) 

+ (Z”(0)^ + t ||M^||g5/2Z”(T)dT. (4.23) 

Jo 2,1 

Now, we define m G Z by 


m 


= inf < p G Z| ^ 2 ^ 7 / 2 11 11^2 < Cob 

[ ^>P J 


for some sufficiently small positive constant Cq. 

Noticing that div(w” + u”) = 0, we get by applying Proposition l3.2l to the first equation of (14.11) 
that 


I‘^”IIlP(B3(2) < I|Mo|Ib 3/2 exp |c (^||Mg||gl(g5/2) + l|M”||gl(B5/2)^ I 


(4.24) 
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and 


I^”IIl“(L2) < II^oIIlZ/ lk”llL”(L“) < ||«o||l“/ 


(4.25) 


from which, we can get 


< 2®'"||flo||t2 fl+ ||«o||i~ + IIM 0 III/ 2 ) +22'"||flo||i» = Nm. 


2,1 


(4.26) 


Inserting (I4.22D (I4.23D into (I4.19D yields that 


\U 11 y 00 / T3l/2'j H~ ||W 




5/2) + I|VP"||^1(b1/2) 


-^(II“”IIl 1(B2(^) ^io II“”IIb 5/2 + ll‘^LllB2/2)Z"(T)iiT 

+ C(Z"(f))^ + C(l+ ||fl"||j„(g3/2p(l+ ||Wo|IbV 2 + ll‘^o|lB3/2)(||ML||gl(B5/2) + ||<|Ib1(b2/2)) 

+ a2«-||s.«'^|lt»(,2)(i + ||«1lt.(g3/2p(z"( 

+ Cf2^"'||S^fl"||g»^g2)(l + ll®”lll~(B3/2p(l + ||Mo|||i/2 + IM 0 III 3 / 2 ) 

+ Cf2^"'||S^fl"||g»(g2)(l + ||«"||£^(b3/2P(||Mo|Ib 1/2 + ||do||g3/2)(Z"(f))'^ 

+ Cf28-||S^««||4„^^,^(l+||««||4^^.3,,p|l^ (4.27) 


Summing the above result up to (I4.16|) . together with (14.31) . (I4.24|) . (I4.25|) and taking e small enough, 
we can finally get 

^ (0 =l|w I|£~(b1_/2) + 11^ IIl“(b5/2) + I|w IIl1(B|(^) 11“^” + 11^-^ 

< j\N„ + ||h"J|j„. + (1 + ||d„||j,,.)||<||j„.)Z"(T)rfT + C(1 + lldoll j,,. + (Z"(f))^)(Z"(())" 

+ (1 + ll«o|lB3/2exp(C||wo||gl/2)e*^ **’*)(! + Hwollgl/s + IMo|lB3/2)(||liL||,,jg5/2j + MlIIl1(b)22)) 

+ f2»"||«„||t.(l + ||<loll*,., 2 exp(C||«„|| + IldolljM + Z"(())(Z"(())* 

+ f2»”||«„||t.(l + ||B„||*„exp(C||«„|| + ||«„|||.„ + ||4I||,,2) 

+ f2»“|Moll‘B2(l + ||«oll‘„ 2 exp(C|| 2 ,„||B.,.)e'z"l')))||rf„||*,.(Z»(f)), (4.28) 

" 2,1 2,1 “2,1 

under the assumption 


I(7 S^ci ||gc«^g3/2^ ^ 2cnI 7 . 
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Gronwall inequality help us to get 


Z”(f) <C3 exp(C3(fNm + ||Mo||gi /2 + (1 + ||rfo|lB3/2 ) |Mo|Ib3 /2)) 

X |c3(l+||do|lB3/2 + Z"(t))(Z"(t))2 

+ C3(l + ||flo|lB3/2exp(C3||uo||gi/2)e’'^ + ||mo|Ib1/2 + ||do|lB3/2) 

l>2i ^->2,1 ^2,1 ^2,1 

\q€Z ) 

+ C3t2®”’||flo|lt2(l + lko|||3/2exp(C3||wo|lBi/2)e^^ ^*^^)(lk'o|lBi/2 + ||rfo|lB3/2 + Z"{t)){Z’\t)Y 

+ C3t2®"'||flo|lt2(l + ll«o||t3/2exp(C3||i/o||Bi/2)e^^"(*)))(l + ||mo|||i/ 2 + ||<^o||j3/2) 

^2,1 2,1 122,1 ^2,1 


+ C3t2®"'||flo|lt2(l+ ll«o||j3/2exp(C3||Mo||Bi/2)e(^"(*)))||(io|lt3/2(Z”(t)) 


D-5/^ 

^2,1 


(4.29) 


where we assume the constant C 3 > C 2 > 1 . 

Applying (I3.5II to the first equation of (I4.1D together with (I4.2I) . (I4.3I) , (I4.24D , (I4.25D and 1 + xe^ > 
for Vx > 0, we deduce that 


\a" - S,„fl'’||£„(g3/2) < ||fl” - S^fl”||fo„^„3/2^ 


iLr(B|f) 

< ^ 23'?/2||A^flo||g2 + ||flo|lB3/2 (expC(Z'’(f) + ||wEIIl1{b5/2)) “ 

q>m 

<Cob+ ||flo||B3/2(Z"(f) + ||Mg||gi^g5/2p exp (^C(Z"(f) + ||wlIIl1(b®(2))) 

<cofo + C2||flo|lB3/2exp(C2||wo||Bi/2)e('^2Z"(d) 2'?/^(l - HA^moIIlz + 2"(f) 

Using the fact that 1 + x < e^, x > 0 and taking 0 < Ti (cq) < (CsNm)^^ so small that 

eC3(l + ||flo|lB3/2)exp(C3) 


(4.30) 


2,1 


X ( +2^''^^IIVoIIl2)] <£o. 


(4.31) 


for some sufficiently small positive constant eo- 


20 






Then it follows from (I4.29I) that 

Z”(f) <eoe^^^ + eCsexp(C3(||Mo|lBi/2 + (1 + ||<^o|Ib3/2)||c1o|Ib3/2)) 

X |c 3 (l + ||<io||j„ + Z”(())(Z”(())" 

+ C3t2»”||«oiry(l + lkoll*„exp(C3||uoll4.x)«'''"''’>)(l!><ol!4.x + MolU™ + Z”(f)) (Z”(f))' 


+ C3f2®'"||flo|lt2(l + ll«o||l3/2 exp(C3||uo|lBi/2)e^^"^‘^^)(l+ ||mo||5i/2 + ||‘^o||53/2) 

" 2,1 2,1 122,1 " 2,1 

+ Cst2^^\\ao\\Ul + \\ao\\ly,exp{Cs\\uoy^^^^^ (4.32) 

122,1 2,1 122,1 I 

provided that 

C2||flo|lB3/2exp(C2||wo||Bi/2)e^‘“2^"^^^^ ( ^ 2'?/^(l - HA^moIIlz + Z”(f) j < Cob. (4.33) 

\qeZ, / 

Letr"(eo) = sup{t € [0, ri(eo)] |Z”(f) < 2eeo}. Without loss of generality, we may assume that £0 

is so small such that 2 eC 3 eo < 1 and 2 e^C 2 £o||flo||B 3/2 6 ^p(C 2 ||wo||Ri/ 2 ) < Cofc.Thenfor t < ri’(eo), we 

^2,1 ^2,1 

get from (14.311) that (14.3311 holds and then 


Z”(f) <££0 + eC3exp(C3(||wo||Bi/2 + (1 + \\do\\g3/2)\\do\\02)){C3{l + UdollB^/z + 2e£o)2e£oZ"(t)) 

|gi/2 + (1+ ||do|lB3/2)||do|lB3/: 


tH1/2 I V J- I II 11 d-5/ ^ J 11 11 R-^/ ^ J J \. '""'A V I II ‘^U 11 ^ 

^2,1 ^2,1 ^2,1 ^2,1 

+ £C3exp(C3(||i/o||gi/2 + (1 + ||iio|lB3/2)ll^o|lB3/2))(^“Ai + fZl2 + tZl3), (4.34) 


with 


Zli = C32^”'||flo||£2(l + e||flo||t3/2 exp(C3||Mo|lBi/2))(ll^o|lBi/2 + II‘^o|Ib 3/2 + 2e£o)(2e£o)^, (4.35) 


"^2,1 


2,1 


"^2,1 


A 2 = C32®'" ||flo 11^2 (1 + e||«o 11^3/2 exp(C3||wo|lBi/2))(l + ||wo||^i/2 + IMo||^3/2 


I UJ/- 
^2,1 


Ri/Z 

^2,1 


^2,1 


and 


(4.36) 


ZI3 = C32^"'||flo 1112 (l + e||«o 11^3/2 exp(C3||Mo||Bi/2))||do||^3/2(2e£o). 

122,1 2,1 122,1 


Taking £0 so small that 


2e C3£o( 1 + ||do|lB3/2 + 2e£o) exp(C3(||Mo|lBJZ + (1 + II^o|Ib|/ 2 ) II^o|Ib3/2)) <1/15 


^2,1 


' 2,1 


and 


Tiieo) = min (Ti(£o), 11,12,13) 


(4.37) 


(4.38) 


(4.39) 
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where 


h < eeo/exp(C3(||Mo|lBi/2 + (1 + ll^o|lB3/2)||rfo|lB3/2))^, (4.40) 

ti < eeo/(^eClA2exp{C3{\\uo\\j^u2 + (1 + IMo|lB3/2)||rfo|lB3/2))) (4.41) 

and 

t3 < eeo/(eCf^3exp(C3(||i/o||Bi/2 + (1 + ||t?o|lB3/2)||do|lB3/2))). (4.42) 

Therefore, we have 

Z"(f) < 4e£o/3, for V < T 2 (eo). (4.43) 

This, together with (14.21) and (14.31) . ensures 

{a",u",d",VP"}neK is uniformly bounded in (4.44) 


STEP 3: (Convergence). The convergence of (w”, d") to {ui, di) readily stems from the definition 
of Besov spaces. As for the convergence of (fl”, u",d"), it relies upon Ascoli's theorem compactness 
properties of the consequence, which are obtained by considering the time derivative of the solution, 
we omit the details here .Moreover, there holds 

I|fl” - Smfl”||B~(B3/2j < 2cob (4.45) 

with the constant Cq being sufficiently small. 

4.2. Uniqueness of Theorem 1.1 

Let {af u\ d\ VP') (with i = 1,2) be two solutions of the system (11.21) , which satisfy (14.451) and 

{a\ u\d\ VP') GCb([ 0 , T]; ^ qQq, T]; n L\[0, T]; ) 

X Cfc([0, T]; B^f) n pi ([0, T]; x ([0, T]; B^^^). ( 4 . 45 ) 

We define 

{Sa, Su, 5d, VSP) = {a^ - - df, VP^ - VP^ 
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so that {5a,5u,5d,V5P) solves 


' dt5a + ■ 'S/da = —5u ■ Vfl^, 

dtSu + ■ V5u + (1 + a^) {VdP - A3u) = F, 

< dtSd - Add + ■ V3d = 3u ■ + Vd^ ■ V3d ■ d^ + V3d ■ Vd} ■ d^ + Vd^ ■ Vd^ • 3d, (4.47) 

div^M = 0, 

, {3a,3u,3d)\t=o = (0,0,0), 


with 


F = - (Jm • Vm^ + 3a{Au^ - VP^) + (1 + a^){Vd} • A5d) + (1 + a^){V3d ■ Ad^) + 3aVd} ■ Ad^. 
Applying Proposition l3.1 1 to the first equation in (I4.47D yields 



(4.48) 


For the momentum equation of (14.11) . we can follow the proof of y] (up to a slight modification) 
to get 



where we have used the following key estimates which are different to Navier-Stokes equations: 



(4.50) 
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||{l+n')(V<i'-AM)||m„_., <C J (l + |lfl'l|B3,0l|V<i'||,„||Airf||„-.dT 


<C 


j(|(l + ll<l'llB3,3)||Vd'||,„||<S<i||„.3iT 

<C /''(1+ ll<l'llBM)l|V«i'||B3,3||id||=« ||M||‘«dT 

JO 2,1 2,1 D2,1 “2,1 

<£||^5d||^l(g5/2) + (1 + ||fl^||g3/2)(l + ||‘^^|lB5/2)||^rf||Bl/2(iT 

<£||l5d||gl(g5/2) + (1 + ||fl^|lB3/2)(l + ||c?^||B3/2||rf^||B7/2)||W||gl/2dT 

(4.51) 


<£pd||gi(g5/2) +C f {l+\\a^\\l,;2)il+\\d^\\l5/2)\M\Bl{2dT 


<£pd||gl(g5/2) + C l^{l+\\a ||g3/2)(l + llBfj2|M llB2/2)||^,d||gl/2dT, 

(4.52) 


\\3aVd^ ■ <cj^ \\3a\\^^n\\Vd^\\L2\\M^\\L-dT 

||dl||g3/2||d2||g3/2rfT. 


(4.53) 


Similar to the estimate (14.51) . we can get from the third equation of (14.471) that 

< x: • V, A,]iid||Bi(B 2 ) + \\3u • Vdi||gi(Bi/ 2 ) 

+ II • V^d • d'^\\jj^j^i/2-^ + II Vdd • Vd^ • d^||gi^gi/2^ + || Vd^ • Vd^ • dd||gij'gi/2^. (4.54) 
Applying Lemma IZ5l 12.71 gives that 

■ V,A,]dd||gi(g2) < ll^rf|lBi/2||M'||B3/2dT, (4.55) 


11 du • Vdi 11 (gl/2) < C ^ 11 du 11 gl/2 11 di 11 g5/2dT 


.7/2 


dr 


<C [ ||‘Jw||g-l2||^w||i//\||di||]^02||rf^ir^^ 

Jo “ 2,1 “ 2,1 “ 2,1 


<£||du||gi(g3/2) +C ||dw||B-l/2||d^||g3/2||d^||g7/2dT, 


(4.56) 
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dr 


dr 


||Vd -VSd-d llj^l^gl/ 2 ^ <C^ ||W||g 3/2 lid ||g 3 / 2 ||d ||g 5/2 

<e 11 dd 11 gi (g5/2) + C ^ 11 dd 11 gl/2 11 d^ 11 3 3,2 11 d" 11 g7/2dT, 

IlVdd • Vdl ■ d2||gl(gl/2) <cj^ ||dd||g3/2||d2||g3/2||d3||g5/2dT 


(4.57) 


< C / JI dd 113 ^2 11 dd 11 ^5^2 11 d^ 11 g3/2 11 d^ 11 3 / 3^2 11 11 g7/2 


|l/2 11^2 


l||l/2 ||jl||l/2 


' 2,1 


'B, 


’ 2,1 


<e||dd||gi(g5/2) + ¥d\\j^in\\d^\\j^3n\\d^\\jj3n\\d^ 

|Vd3 ■ Vd^ -ddll 

L]<cljm ||d^||g5/2dT 

<C f\\M\\ gl/2 lld^llgs/zlld^ll g7/2 dr. 


(4.58) 


(4.59) 


Choosing e small enough and inserting the above estimates into (14.541) . we have 




< 


[ + II^‘^IIb1(2)(I|w^IIb 5/2 + ll^^llB3/2||d^||g7/2 + (1 + ||d^||g3/2)||d3|| 3/2||d3|| 7/2)dT. 

J 0 2,1 2,1 2,1 ^2,1 2,1 2,1 2,1 2,1 

(4.60) 


Combining (I4.49D with (I4.60D , applying Gronwall's inequality and using (14.4611 implies Sa = 3u = 
Sd = 0 for all t G [0, T]. 

5. Global well-posedness of Theorem 1.1 

5.1. The estimate of the transport equation 

In this section, we shall investigate the following transport equation: 

dta + u ■ Va = 0, 


divu = 0, 

fl|f=o = «o- 
l/2/'Tp3\3 n r 1 /'n5/2/ 


(5.1) 


Proposition 5.1. Let u = (w^M3) G L^(B3G(]p.3)) p with divu = 0 and ao G B^(^(R3). 

Then 05.11) has a unique solution a G C([0, T]; so l^at 


l'^llL~(B3f) ^ II«o|Ib 3/2 + C||fl||L~(B3/2)(ll“^llgl^g5/2)l|w''ll^ 


,h\\l/2 


h IIL 




(5.2) 


for any t G (0, T]. 
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Proof : The existence and uniqueness of solutions to dS.lD essentially follow from the estimate 
(I5.2I) for some appropriate solutions to (I5.1|) . For simplicity, we just present the estimate (I5.2|) for 
smooth enough solutions of dS.ll) . Thanks to Bony's decomposition d2.1l) , we obtain 


M • Vfl = T„ • Vfl + TZ{u, Vfl) 


(5.3) 


Applying Ay to d5.ll) and taking inner product with Aja, we obtain 
1 d 

- —||Ayfl ||^2 + {Aj{Tu ■ Vfl),Ayfl) + (Ay(7^(l/, Vfl)),Ayfl). (5.4) 

Using a standard commutator's argument and the energy estimate in 1^, we have 

||Ayfl||L2 <||Ayflo||L2 + C / ( ( || [Ay, Sy/_i] Ay/Vfl||i^2 + || (Sy/_iU - Sy_iw) AyAy/Vfl || £2 )dT 

•^0 \H'\<5 


+ c [\ x: {\\Ajn{u,\/a)\y))dT. 


(5.5) 


I <5 

We first get by applying the classical estimate on commutators and d2.6|) with m = r = oo that 

E (ll[^•/Sy/_l]Ay/Vfl||£l(£ 2 ) + ||(Sy/_iU - Sy_iM)AyAy/Vfl||£i(£ 2 )) 

E (IIS/'-lVl/^llLlyLoo) + ||Sy/_iVl/^||£l(£ooy)||Ay/fl||£c<,y£ 2 ) 

\hr\<5 

+ E (I|S;'-iVm^ - Sy_iVM^||£i(£ooy + ||||Sy/_iVu^ - Sy_i Vw^ || Ay/fl||£c<,y£ 2 ) 

l;-y'l<5 

<C E E (II +2MIV^'^llLi(L~))ll^y'^llL“(L2) 

|y-y'l<5y"<y'-2 

-Cdji ^^^^II®IIl“(b3/2)(II«^II£1^b5^ + 

Thanks to Lemma l2Al and d2.6l) with m = 2,r = oo, we obtain 


|Ay72.(M, Vfl) |l£i(L2) 


E (l|S;'+2V£fl||£oo(£c«,) ||Ay/M^’||£»(£2) + ||Sy/+293fl||L“(L~(L2)) 

j'>j-No 


<C E E 


7 "“llL”(L 2 )l|Ay'W ||l~(L 2 ) 


f>j~No j"<f+l 


+ 2(2/') II A, 


7"«IIl~(L2)| 


<Cdy2-3;72||^| 


L~(b 3/2)U|W ll£l(B5/2yl|- IIb1(bV2) 


3||l/2 


\u 


h\\l/2 


lLi(B|f)/ 


LriLliL?))) 

^y'“^llLi(L2(L”))) 


Taking above estimates into d5.5l) and taking summation for j G Z, we conclude the proof of d5.2l) . 
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5.2. The estimate of d 

Applying Ay to dl.lll , a standard commutator process gives 

dfAjd — AAjd = —Aj{tuVd) — Aj{'lZ{u,Vd)) + Aj{\Vdfd). 

Taking inner product with Ajd to the above equation, we obtain 
1 d 

--\\Ajd\\l 2 +C 2 ^’\\Ajd\\l 2 <C(Ay(f„Vd),Ayd) + C(Ay(7t(M,Vd)),Ayd) 

+ C(Ay(|Vd|2d),Ayd), 


(5.6) 


(5.7) 


where we have used the following fact: 
there exists a positive constant C so that 


— / AAjdAjddx > C2^''||Ayd||^2- 


/R3 


Using a standard commutator's argument and the energy estimate in [SO, we have 

||Ayd|A2 + C2^^'||Ayd|Ai(L2) 

<||Aydo||L 2 + C||^;(|Vd|^d)||£i(i2) + C||Ay(d X Aci)|Ai(L2) 

+ ‘^/( E (||[Ay,Sy/_iI/]Ay/Vd|A2 + ||(Sy/_iM - Sy_iW)AyAy/VdlA2)dT 

•^0 \j-j'\<5 

+ c[\ x: (||A^- 7 t(M,Vd)|A 2 ))dT. ( 5 . 8 ) 

We first get by applying the classical estimate on commutators and (I2.6D with q = r = co that 

E (ll[^h‘S/'-lM]Ay/Vd|Ai(^2) + ||(Sy/_iM - Sy_iw)AyAy/Vd|Ai(£2)) 

<C Y. (l|S;'-lVM'|Ay(,.) 

I/-/'I <5 

+ C ^ (||Sy/_iVr/^ - Sy_iVu^||^i(£c«,) + ||Sy/_iVw^ - Sy_i | Aiy^co)) || Ay/d||^oo(^2) 

l/-/'l<5 

I/-;'I <5 -2 


<Cd 


y2 + ll“^llLi(B3f))- 
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Thanks to Lemma 12.51 and (12.611 with cj = 2,r = oo, we obtain 


\Aj'R{u, Vd) ||l1(l2) 


M|Sy/+2V;,d||£~(L»)||Ay/U ||l“(L2) 

i'>i-No 

+ \\Sf+2hd\\L^(L^(Ll))\\^j'U^\\Lf’{Ll{L^)) 
j’>i~Noj”<j'+l 

+ 2^^' 11 I L“ (£2) 11 A;' 11 £l (L2(L”; 

ni,,3iii/2 ii,,;!|| 1/2 , n,,;; 




and 


lA 




(5.9) 


Substituting above estimates into (15.51) and taking summation for j G Z, we can get 


\\d\\ 


^nw) 


,3/2n + C 


'LbBii) 


7/2n <\\do ||g3/2 + C 


2,1 


I J~(B3/2) 11^11^1(87/2) 


+ C||d|| 


£» (B3/2) ( 11 W 11 ^^5/2^ 11 M 11 + 11 M 


|,3(,5/2)). □ (5.10) 


5.3. The estimate of the pressure 

As is well known, the main difficulty in the study of the well-posedness of incompressible MHD 
system is to the derive the estimate for the pressure term. We first get by taking div to the second 
equation of (11.21) that 


—AP =div(flVP) + div(M • Vm) — div(flAM) — div[(l + a){Vd ■ Ad)] 

=div(flVP) + div;,div;,(M^ (8) u^) + 233div;,(M^w^') — 233(M^div;,M^) — div;,(flAw^) 

— ds{aAu^) — div[(l + fl)V • (Vd O Vd)]. (5.11) 


In what follows, we will give the estimate of the pressure which will be used in the estimates of 

Proposition 5.2. Let a G L“(b3/2(R3)) u = {u^,u^) G n L\{bI^^{M.^)), d G 

Lj(B2i^( 1R^)) n Lj(B2i^(1R^))- Then (15.111) has a unique solution VP G Lj(B 2 3 ^(]R^)) which decays 
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to zero when |x| —> oo so that for all t G [0, T], there holds 

||VP||^i{gi/ 2 ) <CY{a,u,d), (5.12) 

provided that C||fl||£oo|.g3/2^ < 1/2, where Y{a,u,d) is defined by 

Y{a,u,d) =C^||i/ ||£»(bi/ 2 )||w |Il1(b5/2) + II^IIl^(b3/2^(||w |Il 1 (b®//) + 11“ IlipB^f)) 

+ f + (1 + ||“||£~(B3/2j)||c?||g»^g3/2)||rf||Bl(g7/2)^ . 

Proof of Proposition l5.2[ As both the existence and uniqueness parts of Proposition 15.21 basi- 
cally follows from the uniform estimate (I5.12D for appropriate approximate solutions of (15.Ill) . For 
simplicity, we just prove (I5.12D for smooth enough solutions of (15.111) . Firstly, applying Ay to (15.111) 
and using Lemma we have 

+ ||Ay(MMiv,,M'^)||gi(g2) + ||Ay(flAM'A||Ll(L2) 

+ ||Ay(flAM^)||giyg2) + |lAy[(l + fl)V ■ (Vd0 Vd)]||Bi(B2). 

Applying Lemma l2Al and Lemma [2.101 gives 


(5.13) 


||Ay(flVP) ||BiyB 2 ) < dy 2 ^ p ^ II“IIl“(b 3 (^) II^-^IIl|(b 1 ( 2 )' 

2lAy(u^ 0 M^)||b1(B2) < 

2lAy(l/^/^)||Bl(B2y + ||Ay(u3div,,M'AllLl(L2) ^ djl^-P^'^F{u^, 


(5.14) 


(5.15) 


(5.16) 


A/(«AM^)||l 1(B2) + ||Ay(flAM^)||Bl(B2) < djl^ ^^^^II“IIl“(b3//)(||“^||£1(b5/2) + 11“^ |Ib1(b5/2)), (5.17) 


||Ay[(l+fl)V • (Vd 0 Vd)]||Bl(B2) <dy2( P^\l+ ||«|lByo(B3/2p||d||g^yg3/2y||d||gl(g7/2y. 

Taking above estimates (15.141) - (I5.18|) into (I5.13|) . we have 


(5.18) 


II^-^IIl}(BA^) —^|^ll‘^llL“(B3^)l|VP||glygl/2y + ll“^llL~(B^//y||M^'||giyg5/2y +P(w^,W^) 

+ ll“llL”(B|f)(ll“^llLi(B®f) “*■ 11“ llL)(B®f)) + (4 + ll“llL“(B®f))ll^llL“(B|f)ll‘^llLi(B^_(^)^ • 

(5.19) 

So provided that C|lfl||goo(g3/2y < 1/2, we can conclude the proof of (I5.12D . 
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5.4. The estimates of the horizontal component. 


Thanks to the second equation of (I1.2D again, we have 

dtu^ + u ■ - (1 + fl)(Au^ - VhP) = - (1 + fl) V • (Vd 0 V;,d). (5.20) 

Applying Aj to (15.201) and taking inner product with AjU^\ we obtain 

^I|A;W^||l 2 + Ci2^'^||AyM^||^2 

<\\AjuUL^ + I|A;(m • Vu’^)\y + ||Ay((l + fl)V • (Vd O V,d))|A2 
+ ||A^-(flAM^)||^2 + ||A^-((l + fl)V,P)|A2. (5.21) 

Integrating the resulting inequality over [0,t]and using the fact divw = 0, one shows that 

l|A/M^'||Lf’(L2) + 2^^'(l|A;W^||£l(L2) 

~I|A/M()||l2 + II A; (m • Vm'AIIl|(l 2) + II A; ((1 + fl)V • (Vd 0 V;,d)) ||^i(^2) 

+ l|A/((l + «)V;,P)||£i(^2) + ||Ay(flAw^)||^i(L2) 

~I|A/Mo||l2 +2-''||w^ 0 m^IIl1(l2) +2^'||w^m^||l1(l 2) + IIA; ((1 + fl)V • (Vd 0 V;,d)) ||l1(£2) 

+ l|A/((l + fl)V;,P)||Li(L2) + ||Ay(flAM^)||^i(^2). (5.22) 

Therefore, we can get from Proposition l5.2l and the estimates (I5.15ll - (I5.18I) that 

11^^ WipiBif) IIlkb^^^) — II“oIIbi /2 + y(fl,M,d). ( 5 . 23 ) 

provided that C MlpiBlf) < 1/2- 
5.5. The estimates of the vertical component 

In this section, we will use the fact that the velocity field equation on the vertical component is 
a linear equation with coefficients depending only on the horizontal components and a. Thanks to 
fll.lD ?. and divw = 0, we only consider the vertical component of the system to get 

dfU^ + ii ■ Vw^ “ (1 + a){Au^ — d^P) = — (1 + fl)V • (Vd 0 d^d). (5.24) 
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Applying the operator Ay to above equations and taking inner product of the resulting equation 
with AjU^ , we get by a similar derivation of (15.22b that 


~II'^;WoIIl2 + w^ii^ • '^W^)IIl 1(L2) + ll^;((l + ^)33b’)||Li(L2) 

+ ll'^/((l + ■ (Vd 0 93 ^^))||l1(l2) + (5.25) 


Applying Lemma IZSl and Young's inequality gives 


||Ay(u.Vu3)||,j(,.) <C21 A^-(uV)||,1 (,. 


lL2(B3f)l0 llL2(B3f) 


<Cdj2^ ^(||M^||£^|.gl/2) + ||u^||^l(-g5/2p(||M ||£co(gl/2^ 


Li(B|f)^' 


(5.26) 


11 Ay (fl Am )||£1(l2) <Cdj2^ ^ II^IIl“(b3{2) 11“ HlKb^^^)' 

|Ay((l + fl)a 3 P)||glyg 2 ) <Cdy2(“^’/^)||fl||g„(g3/2yY(fl,M,fo) 


(5.27) 


and 


||Ay((l + fl)V • (Vd0 93d))||^i(^2) 

<Cdj2^ + II“IIl“(b3(^))(II^IIl“b®(^)IMIIlJ(b^(2)- (5.28) 

Inserting estimates (15.26b - (15.28b into (15.25b , we have 

11 “ llL“(B 2 f) “*“ 11 “ W L] —W^oW bI{^ 

+ C(||M3||g^(gl/2y + l|M^llLl(B5/2))(||M^||g„(gl/2y + 11 m'" 11 £l (b5/2) ) . (5.29) 

provided that C I1 «|Ib»(b3/2) < 1/2. 

5.6. The proof of Theorem 1.1 

The goal of this section is to present the proof of Theorem 1.1. In fact, given mq ^ 

“0 e B 2 ,i^(R^),do e B2 ,i^(R 3) with ||flo|| g3/2 being sufficiently small, it follows by a similar argument 
in section 4 that there exists a positive time T so that (1.2) has a local solution (a, m, d) with 


a e C([0, T]; ^3/2y 

U e C([0, T]; n jy gl/2) n p], g5/2)^ 

d e C([0, T]; B^f) n p], g3/2) ^ piQo^ p], g7/2)_ (530) 
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We denote T* to be the largest possible time such that there holds (4.1). Then, the proof of Theorem 
1.1 is reduced to show that T* = oo under the assumptions of (11.2b . (11.4b . Let e, rj be small enough 
positive constants, we define T** by 

T** = sup |f G [0, T*) : \\a\\i^^g3^/2^ < 4||flo|lB3/2, 


\\d\\ 


r 3/2^ + l|M''|lLi(Bff) + ll^llLi(Bff) - 4 (II“oIIb1_/2 


^2,1 


11“ ll‘^llL“(B^(h 11“ IlLhBz^/i") II“IIl1(B2";i") 

In what follows, we will prove that j** = j’* under the assumptions of (11.4b , dl.Sb . 
If not, we assume that T** < T*, and for Vf < T**, we get from (15.2b that 


+ \\dyn)+rj^Ho}. (5.31) 


By taking 


we have 


l“llL“(B3f) ^ II“o|Ib 3/2 + C||fl||j„(g3/2)(Ho + (||wg||gi/2 +e)^^^H(5/^). 


C(Ho + (||M^||gi/2 +c) 1 / 2 h 1 / 2 ) < 1 / 2 , 


fl||l»(B3/2) < 2||flo||R3/2 Vf < T**. 


^2,1 


Thanks to (15.29b and (15.31b . we have 

ll“^llL“(Bif) + II“^IIl1(b5(2) <II“oIIb 1// +C||flo|lB3//Wo + C||flo|lB3//l|M^llLi(B^^^^^^ 

+ C(1 + ||«o|lB3//)-f^O + C(II“oIIb1 // + 

where the following estimate has been used: 


(5.32) 


(5.33) 


(5.34) 


(5.35) 


^ ^ II "B!(B5f)l' "L]{Byy 

+ (II“^IIl1(b5/i") + ll“^llL“(Bif))^^^(ll“^llLi(B5/2) + ll“^llLr(B)f))^^^ 

<C(||M3||gU2 + e)i/2H3/2. (5.36) 

One hand, for Vf < T**, while taking 

C\\(Io\\^3/2Ho < 1/8, C||flo|lB3/2 ^ 1/8, C(1 + ||flo|lB3/2)llo — 1/8/ 

C(||M3||gi//+e)i/2H3/2<i/8. (5.37) 
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We can deduce from (I5.35D that 


On the other hand, for t < T**, we get from (15.101) . (I5.23|) that 


(5.38) 


11 “ llL~(B^(h 11^IILj(82,1^) 

< 11 Wq 11 b 1/2 + 11 do 11 b 3/2 + CHo 11 w'" 11 Ll (Bf/2) + C11 flo 11 b 3/2 11 w'" 11 Bl (b5/2) 

+ C||flo|lB3/2(lk*ollB^^^ + e) + CHg + C(1 + 11^0|lB3/2)-fdo + C(||moIIbi/2 (5.39) 

1^21 1 ^ 2,1 ^2,1 ^2,1 

Choosing 


CHo(l + Hq) < 1/8, C||flo|lB3/2 < 1/8, C||flo|lB3/2(||wQ|Li/2 + e) < 1/8, 

r>2i £*2,1 ^2,1 

C(l+ |lflo|lB3/3)Ho^ < 1/8, C(||i/3|Li/2 < 1/8. (5.40) 

r>2,i ^2,1 

we get from (15.391) that 


lir(® 2 f) 11 ^ llj-hfiff) ll‘^llir(S 2 f) 11 ‘^11^1(821^) 


<4Ho/3. 


(5.41) 


Combining with (I5.33|) . (I5.37|) and (15.401) . we can get (15.341) . (I5.38|) and (I5.41|) hold if we take C large 
enough in (11.21) . (11.41) . this contradicts with the definition (15.311) . thus we conclude that T** = T*. 
Consequently, we complete the proof of Theorem 1.1 by standard continuation argument. 


References 


[1] H. Abidi. Equation de Navier-Stokes avec densite et viscosite variables dans Tespace critique. 
Rev. Mat. Iheroani, 23(2):537-586, 2007. 


[2] H. Abidi, G. Gui and P. Zhang. On the wellposedness of 3-D inhomogeneous Navier-Stokes 
equations in the critical spaces. Arch. Ration. Mech. Anal., 204(l):189-230, 2012. 

[3] H. Abidi and M. Paicu. Existence globale pour un fluide inhomogene. Ann. Inst. Fourier (Greno¬ 
ble), 57(3):883-917, 2007. 


[4] S.N. Antontsev, A.V. Kazhikhov and V.N. Monakhov. Boundary Value Problems in Mechanics 
of Nonhomogeneous Fluids, Studies in Mathematics and Its Applications, vol. 22, North-Holland 
Publishing Co., Amsterdam, (translated from the Russian)., 1990. 


33 
























[5] H. Bahouri, J. Y. Chemin and R. Danchin. Fourier Analysis and Nonlinear Partial Differential 
Equations. Grundlehren Math. Wiss., vol. 343, Springer-Verlag, Berlin, Heidelberg, 2011. 

[6] J.M. Bony. Calcul symbolique et propagation des singularites pour les equations aux derivees 
partielles non lineaires. Ann. Sci. Ecole Norm. Swp., 14(4):209-246, 1981. 

[7] K. Chang, W. Ding and R. Ye. Finite-time blow-up of the heat flow of harmonic maps from 
surfaces. /. Differential Geom., 36(2):507-515, 1992. 

[8] R. Danchin. Local theory in critical spaces for compresible viscous and heat-conducting gases. 
Comm. Partial Differential Equations, 26(7-8):1183-1233, 2001. 

[9] R. Danchin. Local and global well-posedness results for flows of inhomogeneous viscous flu¬ 
ids. Adv. Differential Equations, 9(3-4):353-386, 2004. 

[10] R. Danchin. Well-posedness in critical spacesfor barotropic viscous fluids with truly not con¬ 
stant density. Comm. Partial Differential Equations, 32{9):1373-1397, 2007. 

[11] R. Danchin and RB. Mucha. A lagrangian approach for the incompressible Navier-Stokes equa¬ 
tions with variable density. Comm. Pure Appl. Math., 65(10):1458-1480, 2012. 

[12] R. Danchin and P.B. Mucha. Incompressible flows with piecewise constant density. Arch. Ration. 
Mech. Anal, 207(3):991-1023, 2013. 

[13] B. Desjardins. Regularity results for two-dimensional flows of multiphase viscous fluids. Arch. 
Rational Mech. Anal, 137(2):135-158, 1997. 

[14] J.L. Ericksen. Hydrostatic theory of liquid crystal. Arch. Ration. Mech. Anal, 9{l):371-378,1962. 

[15] J. Fan, F. Li and G. Nakamurac. Global strong solution to the 2D density-dependent liquid 
crystal flows with vacuum. Nonlinear Anal, 2014. 

[16] G. Gui and R Zhang. Global smooth solutions to the 2-D inhomogeneous Navier-Stokes equa¬ 
tions with variable viscosity. Chin. Ann. Math. Ser., 30 B(5):607-630, 2009. 

[17] G. Gui and R. Zhang. Stability to the global large solutions of 3-D Navier-Stokes equations. 
Adv. Math., 225 B(3):1248-1284, 2010. 


34 


[18] Y. Hao and X. Liu. The existence and blow-up criterion of liquid crystals system in critical 
besov space. Commun. PureAppl. Anal, 13(l):225-236, 2014. 

[19] J. Hineman and C. Wang. Well-posedness of nematic liquid crystal flow in Arch. 

Ration. Mech. Anal, 210(1):177-218, 2013. 

[20] M. Hong. Global existence of solutions of the simplified ericksen-leslie system in dimension 
two. Calc. Var., 40(1-2): 15-36, 2011. 

[21] T. Huang and C. Wang. Blow up criterion for nematic liquid crystal flows. Comm. Partial 
Differential Equations, 37(5):875-884, 2012. 

[22] O.A. Ladyzenskaja and V.A. Solonnikov. The unique solvability of an initial-boundary value 
problem for viscous incompressible inhomogeneous fluids. Zap. Naucn. Sem. Leningrad. Otdel. 
Mat. Inst. Steklov. (LOMI), 52:52-109,218-219, 1975. 

[23] F. Leslie. Theory of flow phenomenain liquid crystals. Adv. Liq. Cryst., 4:1-81,1979. 

[24] X. Li and D. Wang. Global solution to the incompressible flow of liquid crystal. /. Differential 
Equations, 252(l):745-767, 2012. 

[25] F. Lin and J. Lin. Liquid crystal flow in two dimensions. Arch. Ration. Mech. Anal, 197(5):297- 
336, 2010. 

[26] F. Lin and C. Liu. Nonparabolic dissipative systems modeling the flow of liquid crystals. Comm. 
PureAppl. Math., 48(5):501-537, 1995. 

[27] F. Lin and G. Liu. Partial regularities of the nonlinear dissipative systems modeling the flow of 
liquid crystals. Discrete Contin. Dyn. Syst., A2:l-23,1996. 

[28] F. Lin and C. Wang. On the uniqueness of heat flow of harmonic maps and hydrodynamic flow 
of nematic liquid crystals. Chinese Ann. Math., 197(6):921-938, 2010. 

[29] J. Lin and S. Ding. On the well-posedness for the heat flow of harmonic maps and hydrody¬ 
namic flow of nematic liquid crystals in critical spaces. Math. Methods Appl. Sci., 35(2):158-173, 
2012 . 


35 


[30] P.L Lions. Mathematical Topics in Fluid Mechanics, vol. I: Incompressible Models. Oxford Lecture 
Ser. Math. Appl., vol. 3, Oxford University Press, New York,, 1996. 

[31] Q. Liu and J. Zhao. A regularity criterion for the solution of the nematic liquid crystal flows in 
terms of Boo^oo'riorm,. /. Math. Anal. Appl, 407(2):557-566, 2013. 

[32] M. Paicu and P. Zhang. Global solutions to the 3-D incompressible anisotropic Navier-Stokes 
system in the critical spaces. Comm. Math. Phys., 307{3):713-759, 2011. 

[33] M. Paicu and P. Zhang. Global solutions to the 3-D incompressible inhomogeneous Navier- 
Stokes system. /. Funct. Anal, 262(8):3556-3584, 2012. 

[34] J. Simon. Nonhomogeneous viscous incompressible fluids: Existence of velocity, density, and 
pressure. SIAM J. Math. Anal, 21(5):1093-1117, 1990. 

[35] C. Wang. Well-posedness for the heat flow of harmonic maps and the liquid crystal flow with 
rough initial data. Arch. Ration. Mech. Anal, 200(1):1-19, 2011. 

[36] H. Wen and S. Ding. Solutions of incompressible hydrodynamic flow of liquid crystals. Non¬ 
linear Anal RWA, 12(3):1510-1531, 2011. 

[37] C. Zhai and T. Zhang. Global well-posedness to the 3-D incompressible inhomogeneous 
Navier-Stokes equations with a class of large velocity. arXiv:1410.6300vl [math.AP]. 

[38] T. Zhang. Erratum to: Global wellposed problem for the 3-D incompressible anisotropic 
Navier-Stokes equations in an anisotropic space. Comm. Math. Phys., 295(3):877-884, 2010. 


36 


